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We present an introductory discussion of deep-inelastic lepton-proton scattering as a 
means to probe the substructure of the proton. A resume of QCD is given, emphasizing 
' the running of the coupling constant and the DGLAP evolution equations for the parton 

densities. The determination of parton distributions is discussed and their importance 
for predictions of processes at the LHC is emphasized. Going beyond the pure DGLAP 
regime, we briefly discuss the behaviour of parton densities at low x, and the evidence for 
^ ' non-linear absorptive contributions. 

1 Deep inelastic scattering (DIS) introduced 

High energy electron scattering is an ideal probe of the structure of a composite object. For 
instance, consider the scattering of a beam of electrons on a nuclear target of mass Mat. The 
scattering occurs via the exchange of a virtual photon, see Fig. 1. Since it is virtual, the photon 
is not on its mass shell. That is, its 4-momentum q does not satisfy q 2 = 0. On the other hand, 
a real (ingoing or outgoing) particle or system must be on its mass shell. So the invariant mass 
W of the outgoing system in Fig. 1 satisfies 

W 2 = (p N + q) 2 = M 2 N + 2p N -q + q 2 , (1) 
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Figure 1: Electron- nucleus scattering, where pn and q are the 4- momenta of the incoming 
nucleus and virtual photon respectively, and W is the invariant mass of the outgoing hadronic 
system. The lower three diagrams are a schematic illustration of the cross section for electron- 
nucleus scattering, eN — ► eX, plotted as a function of the scaling variable %n = Q 2 /^Vn • Q at 
three different values of Q 2 . In the lowest plot the wavelength A of the virtual photon probe 
is much less that the nuclear radius R^, and the photon probes a constituent proton of the 
nucleus. 
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where M N and pn are the mass and 4- momentum of the nucleus. It follows that q 2 is negative. 
So we define Q 2 = —q 2 . 

The wavelength of the probing photon A ~ 1/Q- Let us follow what happens as we increase 
the electron energy, so that the photon probe has a shorter and shorter wavelength A. We begin 
with A ^> Rn, where Rn is the radius of the nucleus. In this case the photon sees a "point" 
nucleus and we have elastic electron- nucleus scattering with W = Mn- Thus, from (1), 



where v is the energy loss of the electron. The expression in the laboratory frame shows 
immediately that Q 2 = —q 2 is positive. We sketch the corresponding elastic peak at xn = 1 in 
the first of the three plots of Fig. 1. If we increase Q until A ~ Rn then the outgoing system 
may be an excited nuclear state. Now W > M N and x N < 1, as shown in Fig. 1. 

Finally, if A <C Rn, the photon may probe deep within the nucleus. The nucleus is broken 
up. We have deep (Q 2 3> M N ) inelastic (W 2 3> M N ) electron-nucleus scattering. Indeed, 
the electron may scatter off a constituent proton of the nucleus. In terms of a; at, the resulting 
electron-proton elastic scattering peak will occur at 



but will be smeared out due to the Fermi momentum of the proton bound in the nucleus, see 
Fig. 1. M is the proton mass and A is the number of nucleons in the nucleus. The area under 
the Fermi-smeared peak gives the number of protons in the nucleus, and hence the position 
of the peak determines the number of neutrons. The reduction of the eN elastic peak, with 
increasing Q 2 , reflects the small chance of the A — 1 spectator nucleons all happening to be 
moving in the direction of the outgoing struck proton and reforming the original nucleus. 

Let us increase Q 2 even further. Suppose that protons are made up of three point-like 
quarks, then high-energy electron-proton scattering will simply be a replay of electron-nucleus 
scattering one layer of substructure down. We have an analogous sequence of diagrams to those 
shown in Fig. 1, but with Rn replaced by the proton radius R. Also the scattering probabilities 
should now be plotted in terms of 



where p is the 4-momentum of the proton. The continuous curve in Fig. 2 is the analogue of 
the lowest plot in Fig. 1. It shows the elastic eg-scattering peak Fermi-smeared about x = 1/3, 
together with traces of the elastic ep peak at x ~ 1. If there were no further substructure, this 
curve would persist as Q 2 increases. We would have (Bjorken) scaling; the scattering depends 
only on the ratio x = Q 2 /2p ■ q, and not on the two variables, Q 2 and p ■ q, individually, x is 
known as the Bjorken scaling variable. 
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Figure 2: Schematic illustration of electron-proton scattering as a function of the Bjorken 
scaling variable x = Q 2 /2p ■ q. The proton structure function F 2 is defined in the next section. 
The hadrons N* are excited states of the proton. If the proton consisted of just three valence 
point-like quarks the result would be the continuous curve independent of Q 2 . However with 
increased resolution (higher Q 2 ) the photon may probe one of a pair of sea quarks produced 
from a radiated gluon via g — > qq. Indeed, as Q 2 increases, the proton appears to have more 
and more constituents, which all must share its momentum, and so the distribution skews more 
and more towards small x. This trend from the continuous to the dashed curve is characteristic 
of QCD scaling violations. 

In summary, as Q 2 increases, we first have 'nuclear' scaling with a peak at — 1, then 
violations of scaling, following by 'proton' scaling with a peak at x ~ 1, followed by violations, 
and then 'quark' scaling with a peak at x ~ 1/3. If the quarks themselves had substructure then, 
as Q 2 increases even further, we would enter yet again a region of scaling violations followed 
by another onset of scaling. But history does not seem to repeat itself. Scaling violations are 
observed but these reflect the field theory of quarks and gluons (QCD) with coupling a s . The 
photon "sees" the proton made up of the three quarks (called valence quarks) and an arbitrary 
number of qq pairs (made up of sea quarks). The sea quarks originate from gluons, via g —>■ qq, 
themselves radiated from quarks, see the sketch on the right of Fig. 2. Suppose the photon 
probes a quark carrying a fraction £ of the proton's momentum p. Then for essentially massless 
quarks we have 

(£p + q) 2 = m 2 ~ 0, that is £ ~ Q 2 /2p ■ q = x. (5) 

Consequently as Q 2 increases, more and more partons (that is quarks and gluons) become 
evident which have to share the momentum of the parent proton. Each carries a smaller 
fraction £ = x of the momentum, and we get QCD scaling violations (which, as we will see, 
have the form a s Plog(Q 2 / fi 2 )) as indicated by the dashed line in Fig. 2. On hearing this for 
the first time from Wilczek, one of the discoverers of QCD, a famous experimentalist said 
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(a) NC (b) CC 




Figure 3: Neutral- and charged-current DIS mediated by (7, Z) and W exchange respectively. 

"You expect us to measure logarithms? Not in your lifetime, young man". Yet today the 
high precision DIS data from HERA and earlier fixed-target experiments show exactly the 
QCD logarithmic scaling violations predicted. A collection of plots (which show the scaling 
violations) compiled from these deep inelastic ep scattering data can be found in Section 16 of 
the Review of Particle Properties [1]. Introductory discussions of DIS can be found, for example, 
in Refs.[2, 3, 4, 5, 6, 7, 8]. 



2 The DIS observables: the structure functions 

The DIS process, ep —>■ eX, is shown in Fig. 3(a). We talk of the neutral current (NC) DIS 
mediated by 7 and Z exchange. We also have charged-current (CC) DIS mediated by W 
exchange, shown in the second diagram. Recall that by "deep" we mean Q 2 S> M 2 and by 
"inelastic" we mean W 2 = (p + q) 2 3> M 2 . 

The NC cross section is of the form 

da da 2-nyo? 



xs 



dxdy dxdQ 2 Q 4 



x>i L r w i»> ( 6 ) 



where the sum is over j = 7, Z and 7Z representing photon and Z-boson exchange and the 
interference between them; and where 

( G F Ml \ j Q 2 \ 2 

We see the effects of the 7 and Z propagators, and of the QED coupling a and the Fermi 
coupling Gp. Besides x and Q 2 , associated with the hadronic vertex, we have a variable (y or 
s) which depends the energy of the whole ep system 



V = r ~i = ^ , s = (k + p)' ~ ^. (8) 

p-K V-fc/lab.frame Xy 

Both x and y must lie in the range from to 1. The physical interpretation of y is given in 
(18) below. 
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L^ v is the tensor from the leptonic vertex known in terms of k and k', and W^ v is the unknown 
tensor describing the hadronic vertex. Although is unknown it must be constructed from 
the 4- momenta p, q and the metric tensor g^. For unpolarised DIS, there are three tensor 
forms satisfying the requirements of current conservation q^W^ = q u W flu = 0. In this case the 
general form is 



WV = (- g ^ + ^)F 1 (x,Q 2 ) + ^F 2 (x,Q 2 )-ie^f^F 3 (x,Q 2 ), 
\ q J p ■ q 2p ■ q 



(9) 



where P M = p^ — (p ■ q)q^/q 2 - The observable structure functions, Fi(x,Q 2 ), are functions of 
two scalar variables x and Q 2 which can be constructed from p and q. Note that the last term, 
with a qx p type structure, does not conserve parity. Thus F 3 = if Z exchange is negligible. 
If we insert the general form (9) into (6) and use the known forms of U tv , then, after some 
algebraic manipulation, we find 



da 2na 



2 



{Y + F 2 ± Y.xF 3 - y Fl) (10) 



dxdQ 2 xQ 4 
in the M 2 /Q 2 -> limit, where 

Y± = l±{l-y) 2 and F L = F 2 -2xF 1 . (11) 

A similar expression holds for CC DIS (that is eN — > vX or vN — > eX). For both NC and 
CC processes, the — sign for Y_ is taken for an incoming e + or 17, and the + sign is taken for 
an incoming e~ or v. Complete expressions for the lepton and hadron tensors L^ u , W^, the 
structure functions and the cross sections, including those for polarised DIS, can be found in 
Section 16 of the Review of Particle Properties [1]. 

For the moment let us focus on pure 7 exchange, so F 3 = 0. Even then to determine both 
F 2 and Fl as functions of x and Q 2 we need to measure the y dependence. That is we need 
to perform DIS experiments at a range of ep energies 1 . We will see that Fl = F 2 — 2xFi ~ 0, 
so often the experiments either used QCD to calculate F L or assumed that it was zero, and 
presented results for F 2 (x,Q 2 ). Nowadays, particularly with the advent of high y data, the 
data are presented in terms of the so-called reduced cross section 

a rcd (x, Q 2 ) = F 2 (x, Q 2 ) - (y 2 /Y+)F L (x, Q 2 ). (12) 

The HERA machine at DESY collided 30 GeV electrons head-on with 920 GeV protons, 
giving 

s ~ 4E e E p ~ 10 5 GeV 2 . (13) 

Thus Q 2 ~ xys < 10 5 GeV 2 , and x ~ Q 2 /ys > 10~ 4 for Q 2 = 10 GeV 2 . We see from Fig. 4 that 
the (x, Q 2 ) reach of HERA is about two orders of magnitude better than the earlier fixed-target 
DIS data. So given the data how do we interpret it? Let us start with the Quark Parton Model. 



1 Indccd, the final data runs of HERA were made at lower energies specifically to enable direct measurements 
of Fl to be performed. These should be available for DIS2008. 
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Figure 4: Kinematic domains in x and Q 2 probed by fixed-target and collider experiments, 
shown together with the important constraints they make on the various parton distributions. 



3 The Quark Parton Model 

The basic idea of the QPM is that in the DIS process, ep — > eX, the virtual proton interacts 
with one of the quark constituents of the proton, see Fig. 5(a). We view the process from 
a frame in which the proton is moving very fast so that the relativistic time dilation slows 
down the rate with which the quarks interact with each other. Thus the struck quark appears 
essentially free during the short time (about 1/Q) that it interacts with the photon. As a result 
the ep interaction may be written as an incoherent sum (of probabilities) of scattering from 
single free quarks 

da 

dxdQ 2 ^ Jo ^ jgx ^ J \ydxdQ 2 

where is the probability of finding the quark q in the proton carrying a fraction £ of its 

momentum. The electron-quark cross section has the form 

S(x-Q, (15) 



?/.'*^(^)- (14) 



da eq 2ira 2 e 2 q ( s 2 + u 2 " 



dxdQ 2 s 2 V t 2 
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Figure 5: (a) DIS via the Quark Parton Model; the subprocess eq — > eq occurs at cm. energy 
y/I. (b) eq scattering with equal (opposite) helicities occurs with weighting 1 ((1 — y) 2 ) due to 
angular momentum effects. We also show the centre-of-mass scattering angle 9 and note that 
for essentially massless fermions all the particles have 3- momenta of magnitude \k\. 



where s, i and u are the Mandelstam variables for the eq — *> eq subprocess. Two sets of 
alternative expressions are 

s = (xp + k) 2 ~ 2xp ■ k ~ xs, s ~ Ak 2 , 

t = -Q 2 ~ -xys i = -2k 2 (1 - cos 0) (16) 

u — — s — t — —x(l — y)s u = —2k 2 (1 + cos 0), 

where \k\ and 9 are the magnitude of the e, q three-momenta and the scattering angle in the eq 
centre-of-mass frame. If we insert the first set into (15), then (14) becomes 



da 2na 



E f V e\ [l + (l-yf] 5(x-t). (17) 

g JO J 



dxdQ 2 Q A Y 7 ° 

Insight into the y dependence is obtained by comparing the two sets of equations for s, i and 
u. We see that 

y = -(l-cos0), (18) 

so y = corresponds to forward scattering and y — 1 to backward scattering. If e and q have 
opposite helicities then there can be no backward (0 = n) scattering by the conservation of J z . 
This is the origin of the weighting (1 — y) 2 in Fig. 5(b). Crucial to this argument is the fact 
that at high energies (E 3> mf ermion ) the fermion helicity is conserved at a gauge boson vertex. 

If we re- write the QPM formula (17) in the form 



da 27ra 2 



E f l dfif q {£)elxY + 5{x-0, (19) 



dxdQ 2 xQ A g Jo 

and then compare with the general structure function formula (10), assuming only 7-exchange, 
we obtain 

F 2 = 2xF l = I' d U q (0 %e 2 q 5(x - = £ e 2 q xf q (x). (20) 

n J „ 



The first equality (i.e. F L = 0) is known as the Callan-Gross relation, and holds because the 
quarks have spin 1/2. If the quarks had spin 0, then i*\ would have been 0. Also notice that 
in the QPM the structure functions scale, that is have no Q 2 dependence. 

We have noted that the proton is made of valence quarks (uud) and sea quarks in gg-pairs. 
When probed at a scale Q all quark flavours with m q < Q are active. Usually the flavour is 
used as a shortened notation for a parton distribution. So, for example, 



f u (x) = u(x) = u v (x) + u sca (x), 
fu{x) = u(x) = u sea (x). 



(21) 



We therefore have flavour sum rules 



/ (u — u)dx = f u v dx = 2; / (d — d)dx = \ d v dx = 1. (22) 

JO JO JO JO 

The structure function measurements can be used to reveal the quark flavour composition 
of the proton. From (20) we have 

/4 1 , 1 4 1-1 \ 

2 = x {Q u+ Q d+ r + - + 9 u+ 9 d+ r + -)- (23) 

Using isospin invariance it follows that the neutron structure function is 

/4 1 1 4-11 \ 

F* n = x \-d + -u + -s + ... + -d+ -u + -s + ... . (24) 
2 V9 9 9 9 9 9 / V ; 

Similar formulae can be obtained for CC DIS. For the ep — > uX processes we have 

dxdQi - 2irx{Q* + Mj v ) ^ * Y ~ xF * V Fl Y (25) 
Let us consider e~p — > uX (or up — > e + X); here the basic subprocesses are 



e u — ► ud, with "same" helicities 
e~d^b / u, with "opposite" helicities. 



(26) 



If we now recall that Y± = 1 ± (1 — y) 2 , and use the helicity diagrams of Fig. 5(b), then it 
follows that 

F™~ = 2x(u + d + c + s...), xF™~ =2x{u- d + c-s...). (27) 

For e + p — > vX (and up — > e~X) the structure functions are obtained by the flavour interchanges 
d <-> u, s <-> c, while those for the neutron are obtained from those of the proton by the 
interchange u <-> d. Thus at large x, where the valence quark distribution dominates, we have 

a cc (e~p) ~ xu v , and a cc (e + p) ~ (l-y) 2 xd v . (28) 



8 




Figure 6: An historical plot of SLAC NC- and fixed-target muon neutrino CC-DIS data showing 
agreement with the QPM relations (29) and (30), but also indicating that only 50% of the 
proton's momentum is carried by quarks. 



It is informative to compare F 2 (N) = (Fg + F 2 ")/2 for NC and CC DIS. For NC we have 
from (23) and (24) 

F 2 (eN) = —x(u + u + d + d + ...), (29) 

lo 

whereas for CC processes, vN — > /j,X, it follows from (27) that 

F 2 w (isN) = x(u + u + d + d+ ...). (30) 

An experimental comparison of DIS data in the early 1970s is shown in Fig. 6. The good 
agreement with the QPM relations is evident, but the area under the curve 

f 1 F 2 (vN)dx = ( 1 y^xq{x)dx ~ 0.5, (31) 

shows that only 50% of the proton's momentum is carried by quarks. It provided the first 
(indirect) evidence for the existence of the gluonic component of the proton. 

Before the advent of QCD there was a big puzzle. In DIS the struck quark appears to act as 
if it were free inside the proton. Yet it is never seen. No matter how hard it is hit, a free quark 
never emerges. It is confined within the proton. To say that a quark acts as if it were totally 
free, in a deep inelastic scatter, is not quite correct. We need to allow for the interactions 
of quarks and gluons (QCD), and to see how this improves the QPM description of DIS. At 
the same time, we will see that QCD appears to be able to solve the big puzzle. It offers the 
possibility of an explanation of quark confinement. 
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(a) QED 



(b) QCD 




Figure 7: The lowest-order Feynman graphs for (a) the QED interaction between charged 
electrons and (b) the QCD interaction between coloured quarks via a coloured gluon, (RB). 

4 Resume of QCD 

Colour was first introduced to overcome the statistics problem in the quark model of hadron 
spectroscopy. Take for example the hadron A ++ , which was discovered as a resonance in 7i + p 
scattering. It has spin 3/2 and is made of three u quarks, which can be in a state with all their 
spins parallel. But quarks are fermions and, by the exclusion principle, it should not be possible 
to have three identical quarks in the same state. A way to overcome the anomaly was soon 
proposed, which only later was found to have profound implications. The idea is to give quarks 
an additional attribute, colour, which can take three possible values, say red, green and blue. 
Hadrons are postulated to be colourless or, to be precise, colour singlets of the group SU(3) 
constructed from the fundamental colour triplet of quarks (g#, qc, <?b). Baryons (qqq) and 
mesons (qq) are clearly allowed, but single free quarks are forbidden, since they carry colour. 
In effect the puzzle of quark confinement (the experimental absence of free quarks) has been 
replaced by the puzzle of why colour should be confined. 

Recall that the electromagnetic force between charged particles is mediated by the exchange 
of photons, Fig. 7(a). The strength of the quantum electrodynamic (QED) interaction is 
determined by the charge of the particles. QED can be obtained from a remarkably simple 
symmetry principle: invariance of the theory under local phase transformations of the fields for 
the charged particles. Local means that we can arbitrarily vary the phase from space-time point 
to point. Since phase transformations commute with each other and form a U(l) symmetry 
group, we say that QED is a U(l) Abelian local gauge theory. The particles which emerge 
naturally from the theory to compensate for the phase differences from point-to-point or, in 
other words, to ensure the local gauge (i.e. phase) invariance of the theory, have zero mass and 
spin 1. These "carriers" of the electromagnetic force, the so-called gauge bosons, have exactly 
the properties of, and can be identified with, the familiar photon. 

In analogy, in 1972, a local SU(3) gauge theory, quantum chromodynamics (QCD), was 
proposed as the theory of the strong interaction. In QCD the interaction is mediated by 
the exchange of zero-mass, spin-1 gluons between coloured quarks, Fig. 7(b). In QCD there 
are three different colour charges (red, green and blue) which have to be conserved, so the 
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most general phase transformation is slightly more complicated. To be precise, QCD is based 
on invariance under the non-Abelian SU(3) group of local phase transformations among the 
triplet of colour charges, q = (q^, qc, qB)- The gluons themselves carry colour. In fact, eight 
different colour combinations of gluon are required to neutralize all possible phase differences: 
one colour combination is made explicit in Fig. 7(b). 

Note that the gluons have zero-mass and therefore infinite range, and yet the strong force 
between hadrons has such a very short range. Indeed, it is ironic that the nuclear force, where 
it all began, is now relegated to a residual colour interaction, between colour neutral hadrons. 
The binding of colourless protons and neutrons into nuclei is similar to the van der Waals force 
which binds electrically neutral atoms into a molecule. Since colour is confined, the nuclear 
force must be short range and confined to hadronic dimensions. 



5 The running QCD coupling 

The most crucial feature of QCD is the dependence of the QCD coupling, a s = g 2 /4ir, on Q 2 . 
At first sight, it appears that a dimensionless 2 QCD observable R must, for energies Q ^> m q , 
be independent of Q 2 . The only energy scales in the QCD Lagrangian are the quark masses, 
and since the relevant ones are very light we would expect this scaling property to set in at low 
Q 2 . However this argument is not true in a renormalizable field theory like QCD (or QED). A 
scale enters when we use perturbation theory to calculate the observable 

R = £c n <, (32) 

n 

since we encounter (loop) Feynman diagrams which diverge logarithmically. We need to renor- 
malize (reparameterize) the theory, which introduces a renormalisation scale /i. As a conse- 
quence we find that the dimensionless observable R no longer scales, but has logarithmic scaling 
violations, that is, it has the functional dependence R(log(Q 2 / fi 2 ) , a s (fi 2 )). 

First we discuss the QED coupling. Vacuum polarisation effects (i.e. polarised e + e~-pairs) 
screen the bare electron charge. The screening is least at short photon wavelengths, which 
causes the QED coupling, a = e 2 /4iT, to increase with the energy of the photon. The situation 
is shown in Fig. 8, which also shows the relevant Feynman diagrams. Summing up these 
diagrams we obtain, at large Q 2 , 

2 



a(Q 2 ) = a 



(33) 



1 " ^g M2 

where a cut-off, M, on the loop momentum has been introduced to prevent an infinite contri- 
bution. We may eliminate the dependence of a on this arbitrary parameter M by introducing 
the renormalisation scale fi. From (33) we have 

1 1 1 , Q 2 , 1 1 1 , fi 2 

— -^logTT^, and ——r = — - —log— . (34) 



a(Q 2 ) a 3tt to M 2 ' ' a(n 2 ) a 3tt to M 2 ' 



2 The structure function F 2 is such a dimensionless observable. 
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Figure 8: The vacuum polarisation effects which cause the QED coupling, a = e 2 /4ii, to run. 
The shorter the wavelength of the probing photon the more of the bare electron charge it sees. 
QED determines the running of the coupling a, but experiment fixes the normalisation, which 
is traditionally given in terms of a(0) ~ 1/137. 



The M dependence can be eliminated by subtracting these two relations. In this way we obtain 

a(Q2) = h^vIWp) ' (35) 

In effect, the infinities of the theory have been removed at the price of introducing a renor- 
malisation scale fi. QED predicts the "running" of a, but experiment is needed to predict its 
absolute value. 

Note that, due to basic properties of gauge field theories, the ultraviolet divergences of the 
Feynman diagrams of Fig. 9 mutually cancel, via the so-called Ward identities. This is just as 
well, because it ensures that the renormalized charge of the electron, muon,... remain equal. 

Turning now to QCD we have a new vertex to consider, the triple-gluon vertex, which 
arises since the gluons themselves carry colour charge. This changes everything, as Fig. 10 
shows. There is a new vacuum polarisation diagram with a gluon loop, which antiscreens the 
colour charge, which dominates the screening arising from the rif quark-loop diagrams 3 . As a 
consequence the — 1/37T in (35) becomes +b , with 



rif is the number of active quark flavours, that is the number of quarks with m q < Q. 
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Figure 9: The Ward identities ensure that the ultraviolet divergences of these diagrams mutually 
cancel. 
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Figure 10: In addition to the 'screening' quark loop, QCD has 'antiscreening' from the gluon- 
loop diagram, which arises from the non-Abelian nature of the SU(3)-colour gauge group which, 
in turn, allows a triple-gluon vertex (as well as a quartic-gluon coupling, see Fig. 11). As a 
result the running of the QCD coupling constant, a s , is the 'opposite' of QED. It decreases 
with energy, allowing the use of perturbation theory at high energies. 



where 4 

07T lZ7T 

At a stroke, the non-Abelian nature of QCD has solved the puzzling dilemma of the quark 
model. The big puzzle was that when the proton is hit hard in DIS, the quarks act as if they 
are essentially free; and yet no free quark has ever been seen - they are confined within the 
hadron 5 . This asymptotic freedom and infrared slavery is precisely what the running of the 
QCD coupling indicates, see the a s plot in Fig. 10. 

The ultraviolet divergences of the QCD diagrams analogous to those in Fig. 9 cancel due to 
the Slavnov- Taylor identities of QCD. Moreover for a gauge theory the equality of the qqg and 
ggg couplings is preserved by renormalization. 

Let us return to our dimensionless observable R, which due to renormalization becomes the 
function i?(Q 2 //i 2 , a s {ix 2 )). However R cannot depend on the choice of renormalization scale, 



4 The — 1/37T in (35) becomes — 1/67T for each quark loop in (36) due to the convention used to normalize the 
SU(3) matrices. 

5 Confinement still has to be proven. Lattice QCD is the technique to describe physics in the strong coupling 
regime. 
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Figure 11: The diagrams which specify the two- loop /5-function. Note that besides the triple- 
gluon vertex, QCD contains a quartic gluon coupling which gives rise to the final diagram. 



so we have a renormalization group equation (RGE) 

dR _ / d da s d 

dlog/i 2 \ <91og/x 2 (91og/i 2 da s 

It can be shown that the solution of the RGE gives 

R(Q 2 //2 2 } « s (/x 2 )) = R(l, a s (Q 2 )). (39) 

That is the running of a s determines the Q dependence of R. In general, the running is 
expressed in terms of a /^-function, defined by 

dajdhgfi 2 = (5{a s ). (40) 

So far we have introduced the /3-function at one-loop, in which we have summed up the leading 
logs, that is all the (a s \og(Q 2 / fi 2 )) n contributions. From (36) it is easy to show that this gives 

P(a 8 ) = -b a 2 . (41) 

The two-loop /3-function sums up the next-to-leading logs resulting from the two-loop diagrams 
shown in Fig. 11, which are found to give an term with coefficient b\ = (153 — 19nj)/247r 2 . 
The /9-function then becomes 

(3(a s ) = -b a 2 s -h a*. (42) 

Before returning to DIS and the structure of the proton, let us make a few more notes about 
the properties of the QCD coupling a s . First the coupling at two- loops is the solution of the 
transcendental equation 

> L_ + i^m 

a s (Q 2 ) a s (fi 2 ) b \a s {ii 2 ) 

Next, the value of a s depends on the renormalization scheme. Since in two different schemes the 
values are related by a' s = a s (l + ca s ), it follows that b and &i are scheme independent, whereas 
the higher coefficients are not. Nowadays, most calculations in fixed-order QCD perturbation 



R = 0. 



(38) 



b + bia s (/2 2 
b + M S (Q 2 



bolog—. 
fi 2 



(43) 
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Figure 12: The 0(a s ) diagrams which contribute to the proton structure function. On the first 
diagram we show new variables, the longitudinal momentum fraction y carried by the quark 
and the transverse momentum k t of the emitted gluon, which must be integrated over. 



theory are performed in the so-called modified minimal subtraction (MS) scheme. Thirdly, note 
that the coefficients bi depend on the number of active flavours, rif. As Q increases through 
a flavour threshold we will need to ensure the continuity 6 of a s . For example at the 6-quark 
threshold we will require a<;(m 2 ,4) = a s (m 2 ,5). Finally, recall that perturbative QCD tells us 
how the coupling varies with scale, but not the absolute value itself. The latter is obtained 
from experiment. Traditionally the value is quoted at Q = Mz in the MS scheme for nf = 5; 
the current value, determined from many independent experiments is 

a s (M z ) = 0.117 6 ± 0.002. (44) 

Lattice QCD can also predict the value of the coupling. It is encouraging that these very 
different determinations are found to be consistent with each other. 



6 Running parton densities: DGLAP equations 

In terms of QCD perturbation theory, the QPM formula (20) may be regarded as the zeroth- 
order term in the expansion of F 2 as a power series in a s . To include the 0(a s ) QCD corrections, 
we have to calculate the photon-parton subprocess diagrams shown in Fig. 12. The QPM of 
Section 3 is the first diagram in the second |...| 2 . Due to the propagator of the virtual quark, 
the first diagram in the first |...| 2 is proportional to (yp — k)~ 2 oc (2p- fc) -1 , for massless quarks. 
It therefore has a collinear divergence when the gluon of 4-momentum k is emitted parallel to 
the incoming quark of 4-momentum yp. 

The partonic approach is only valid if the second diagram in the first |...| 2 can be neglected. 
Then, the emitted gluon can be considered as part of the proton structure. It turns out that 
both diagrams are required to ensure gauge invariance, but that the second only plays the role 
of cancelling the contributions from the unphysical polarization states of the gluon. Adopting 
a physical gauge, in which we sum only over transverse gluons, only the first diagram remains. 

6 Bcyond two-loops there are discontinuities in a s at the flavour thresholds. Of course, the observables are 
continuous, since the above discontinuities are cancelled by ones occurring in the coefficients functions, see (75) 
below. 
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(45) 



To evaluate its contribution we need to sum over all the possible values of the new variables, 
y and k t , that are introduced to describe the gluon. First we write down the answer and then 
explain its structure 7 . The result is 

x q Jx y 

where P and C are known functions. These will become the universal parton splitting functions 
(here P = P qq describes the q — > qg splitting) and the process dependent coefficient function 
C. The 5-function term in (45) is the zeroth-order QPM contribution with y = x that we 
derived in Section 3. The first order term in a s comes from the first diagram in Fig. 12. A 
straightforward application of the Feynman rules shows that 

= (46) 

and that the log(Q 2 / /i 2 ,) originates from the integration over the gluon (bremsstrahlung) trans- 
verse momentum spectrum 

^dk 2 , (Q^ 




(47) 

/W 

where the upper limit is set by the virtuality of the photon which scatters off a quark of 
transverse size 1/Q. Really the lower limit of integration should be set to zero. We have 
therefore arbitrarily cut-off the integral at some scale /j,q. How do we make sense of (45)? 

Inspection of (45) shows that, after including the 0(a s ) contribution, we may replace f q by 
a well-behaved running parton density 

'.<*• " 2 > = '•<*) + I 7^) I ( p (?) log @) + c ') • < 48 » 



such that 



F 2 (x,Q 2 ) ^f^dy 



x Jx y 



q J x V 




2 



2 



(49) 



where the division of the known function C into C\ + Ci depends on the choice of (factorization) 
scheme. 

The dependence of f q (x,fi 2 ) on the (non-perturbative) scale /xo can be eliminated in an 
analogous way to the dependence of the coupling a on the cut-off M, see (34) and (40). From 
(48) we obtain 

8 -^ 1 = ^f'^My,^p(-), (50) 

dlogfi 2 2?r Jx y \yj 
which describes the evolution of the parton density with fx 2 . This is known as the DGLAP 
evolution equation [9]. In effect we have absorbed all the collinear infrared sensitivity into a 
well defined running parton density f q (x, /i 2 ). We cannot use perturbative QCD to calculate the 
absolute value of f q (x,fi 2 ), but we can, via the DGLAP equation, determine its \i dependence. 

7 This will be a prolonged discussion. Only in Section 8 will we emphasize the important "factorization 
theorem" structure of the formula. 

8 This follows because everything else in (49) is well-behaved. 
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7 Further discussion of the DGLAP evolution equations 



Our 0(a s ) treatment is incomplete. A hint that this is so, is the presence of a soft divergence, 
which arises when the energy of the emitted gluon tends to zero, in addition to the collinear 
divergence. This is reflected in the z = 1 singularity of P(z) = P qq (z), see (46). At this point 
it is crucial to include the virtual gluon diagrams of Fig. 12. The (0(a s )) contribution is the 
interference of these three diagrams with the QPM diagram. This contribution is also singular 
at z = 1. It turns out that the singularity exactly cancels the z — 1 singularity present in the 
real contribution. After the cancellation of the singularity there remains a residual 5(1 — z) 
contribution from the virtual diagrams. Instead of calculating this contribution explicitly, there 
is an easy way to see what it must give. It must be such to satisfy 



Jo 



P qq {z) dz = 0, (51) 



which expresses the fact that the number of valence quarks is conserved during the evolution. 
The virtual diagrams regularize the 1/(1 — z) singularity in P qq so the constraint holds. This 
modification to P qq can be expressed in terms of the so-called "+ prescription" for regularization 
in which 1/(1 — z) is replaced by 1/(1 — z)+ defined so that 

where (1 — z) + — (1 — z) for z < 1. Ensuring the constraint (51) gives 

p « = l0r + + 2S(1 - z) - (53) 

Our 0(a s ) treatment is still not complete. In addition to the ^q — > gq subprocesses shown 
in Fig. 12, at 0(a s ), we need to include the 7*7 — > qq processes. Then the DGLAP evolution 
equation (50) for the quark density q = f q becomes 

= % (P-a, + p.®,) (54) 

where g = f g is the gluon density, and P qq = P is the q — > q(g) splitting function of (53). It 
can be shown that the g — > q splitting function 

P qg = l -{z* + {l-zY). (55) 

In general P a b describes the b — > a parton splitting. Also in (54) we have used <S> to abbreviate 
the convolution integral 

P®f= f-mpi-)- (56) 
J* y \yj 
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Clearly we must also consider the evolution of the gluon density 

dg(x, Q 2 



<91ogQ 2 



= ^ ^E p M®(ft + ft) + p gg®9j, (57) 



where the sum is over the i quark flavours, and where the q — ^ g and g — > g splitting functions 
can be shown to be 

P gq = P qq (l ~z) = - 1 + {1 ~ Z) \ (58) 

p » = 6 (¥ + (T^y: + * - *)) + (y-t) * - 2 »- < 59 > 

Here the coefficient of the 5(1 — z) term can be obtained from the constraint that all of the 
momentum of the proton must be carried by its constituents 



j'dzz^q^Q^+q^Q^+g^Q 2 )^ = l 



(60) 



for all Q 2 . 

It is convenient to introduce flavour singlet (£) and non-singlet (g NS ) quark distributions: 

2 = £(<& + *)■ (61) 

i 

An example of a non-singlet is the up valence distribution 

u v = u — u. (62) 

Non-singlet evolution satisfies (50) and decouples from the singlet and gluon evolution equa- 
tions, which are coupled together as follows 

d f E\ _ a s f P qq 2n f P qg \ ^{ Z\ ^ 



OlogQ 2 V 9 J 27T V Pgq Pgg J \ 9 J ' 

In general the splitting functions can be expressed as a power series in a s 

P ab (a s , z) = P™(z) + a s P^°(z) + a 2 s PT L °(z) + - (64) 

where the NLO expressions were computed in the period 1977-80 and the NNLO in the period 
ending 2004. Leading order (LO) DGLAP evolution, which we have outlined, sums up the lead- 
ing log contributions (a s \ogQ 2 ) n , and next-to-leading order evolution includes the summation 
of the a s (a s \ogQ 2 ) n ^ 1 terms. 

If we are given the x dependence of the parton densities at some input scale Qq then we may 
solve the evolution equations to determine them at higher Q 2 . Frequently this is performed 
simply by step-by-step integration up in Q 2 . 
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An alternative procedure is to rewrite the equations in terms of moments, which for an 
arbitrary function f(z) are defined as 

/<»> = C-z n f{z). (65) 

JO z 

If we now multiply the DGLAP equation (50) by x 11 ^ 1 , and integrate over x, we obtain 

f) rl (y rl rl 

x n - l q^{x,Q 2 )dx = y^OssM^dy J q z n - l P qq {z)dz (66) 

using x = yz. That is, the evolution equation then turns into an ordinary linear differential 
equation for the moments, 

r\ (n) 

ggNg _ ^£p(n) (n) (f]7] 



For fixed a, the solution is 



{n \Q 2 ) = c n exp ( 7 (n) logQ 2 ) = c n [Q 2 y (n \ (68) 



<?NS 



where 7 (n) = a s is known as the 'anomalous dimension'. If we incorporate the running 

of a s , (36), then it is easy to show that 

9ns (<3 2 ) = cn [a s {Q 2 )V {n)/2 ^- (69) 

This is the LO behaviour. In analogy with (38) and (39), the general result may be obtained 
from the RGE 

<V " 1 / ° +/% S )^ + 7 ( " ) (« S )V ) = 0- (70) 



dlog/i 2 y (91og/i 2 da 
which can be shown to have the solution 



qV>{(?lf y *.tf)) = q^(l,a s (Q 2 )) exp f P ( f > ^^-da.) ■ (71) 



In addition to (67) we have 



Once we have the analytic solutions of these equations for the moments, we can obtain the z 
distributions of the partons by the inverse Mellin transforms 

1 fC+iOD 

f t (z,Q 2 ) = — / dnz- n f t (n,Q 2 ), (73) 
where the contour is to the right of all the singularities of the integrand. 
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8 Observables: the factorization theorem 



We return to the equation for F 2 (x, Q 2 ), eq.(49). We have described how the collinear singu- 
larities of the formula have been swept into well-defined running parton densities, fi{y,fi 2 p), 
evaluated at some (factorization) scale 9 /xp in the perturbative region. A convenient choice is 
to set /if = Q, so that the \og(Q 2 /[i F ) term disappears. We then have, including the 75 — > qq 
contribution, 



F 2 (x,Q 2 ) = ^e 2 ! 1 ^- 



2\ " I 



(74) 



where the C 2 ,i are the coefficient functions for the observable F 2 . Although all the collinear 
singularities are absorbed by the running of the f\, recall that the prescription is not unique. 
We can add any finite term. So we must specify a scheme. The MS factorization scheme is 
favoured. It was mentioned at the end of Section 5 as also the choice of renormalization scheme. 

We can generalize this result to describe the structure functions of all DIS processes. For the 
structure functions F a , describing the deep inelastic processes £ + p — * £' + X, the factorization 
formula, which holds to all orders in perturbation theory, has the generic form 

F a (x,Q 2 ) = E f-hiy.Q^CaA-^siQ 2 )) + o(^). (75) 

i= q ,g, g Jo y \y J \ Q J 

The final term denotes non-perturbative contributions, such as hadronization effects, multipar- 
ton interactions etc. For sufficiently high Q 2 these effects are negligible, and the expression for 
the observable factorizes into 

• universal parton densities (of the proton), f iy which absorb the long distance collinear 
singularities. They cannot be calculated in perturbative QCD, but their Q 2 dependence 
is calculable using the DGLAP evolution equations, in which the splitting functions are 
calculable as power series in a s . 

• coefficient functions, C a>i) which describe the short distance subprocess. They are calcu- 
lable from perturbative QCD as a power series in a s , but are unique to the particular 
observable, F a . 

The factorization is displayed visually in Fig. 13 

A similar factorization applies to inclusive 'hard' hadron-hadron collisions. For instance, 
consider the LHC process 

P(pi) + P(p 2 ) - H(Q,..) + X (76) 



9 The subscript F is added to distinguish it from the renormalization scale introduced in Section 5. In practice 
these scales are often chosen to be equal. 
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i= q i= g 




Figure 13: Schematic picture of the factorization theorem for a deep inelastic structure function 
of the proton. 



where H denotes the triggered hard system, such as a weak boson, a pair of jets, a Higgs boson 
etc. The typical hard scale Q could be the invariant mass of H or the transverse momentum 
of a jet. Then according to the factorization theorem the cross section is of the form 



a = 53 / dx 1 dx 2 fi(x 1 ,fJ l 2 F )f j (x 2 ,fJ l 2 F ) 6- ij (xip 1 ,x 2 p2,Q---;fJ?F), 



(77) 



where typically x min > Q 2 / s where s = (pi +p 2 ) 2 - For the production of a system H of invariant 
mass M and rapidity y, the momentum fractions Xi j2 = Me ±y / y/s. The fi and a depend on 
the renormalization scale fiR via a s (/x|.). For instance 



a 



(78) 



m=0 



where LO, NLO... correspond to n—1,2...; note that, for example, k=0,2,.. for W, dijet,... 
production. We should work to the same order in the series expansion of the splitting functions. 
In practical applications it is usual to choose hf = Hr ~ Q and to use variations about this 
value to estimate the uncertainty in the predictions. Of course the physical cross section a 
does not depend on the scales, but the truncation of the perturbative series brings in scale 
dependence. If we truncate at order a™, then the uncertainty is of order a™ +1 . 



9 Global parton analyses 

Two groups (CTEQ [10] and MRST [11]) have used all available deep inelastic and related hard 
scattering data involving incoming protons (and antiprotons) to determine the parton densities, 
fi, of the proton. The procedure is to parametrize the x dependence of fi(x,Ql) at some low, 
yet perturbative, scale Qq. Then to use the DGLAP equations to evolve the fi up in Q 2 , and 
to fit to all the available data (DIS structure functions, Drell-Yan production, Tevatron jet and 
W production...) to determine the values of the input parameters. In principle there are 11 
parton distributions (w, u, d, d, s, s, c, c, b, b, g). However m c , m;, 3> Aqcd- So c = c and b = b 
are calculated from perturbative QCD via g — > QQ. Also the evidence from neutrino-produced 
dimuon data, uN — > n + fi~X, is that 10 s ~ s ~ 0.2(w + d) at Q 2 ~ 1 GeV 2 . 

10 Analysis of NuTeV data for v and v beams indicates some x dependence of the factor "0.2" , and that s > s 
for x - 0.01. 
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A common choice of parametrization of the parton densities is 

xf(x,Ql) = A(l-xfx a (l + e y /x + 'yx) (79) 

with up to five parameters (A, a, (3, e, 7) for each parton. Three of the A's are determined from 
sum rules. The input partons must satisfy the two valence quark sum rules 

f 1 dx(u - u) = 2, C dx{d- d) = 1, (80) 

JO JO 

and also we must satisfy the momentum sum rule (60). 

We can obtain some idea of what to expect for the values of the $ parameters from the 
spectator counting rules. As x — > 1 physical arguments indicate that 

f(x) - (1-x) 2 "-- 1 (81) 

where n s is the minimum number of spectator quarks which share between them the residual, 
vanishingly small momentum of the proton. The greater the number of spectators, the smaller 
the chance of producing a parton with a large fraction of the proton's momentum. For a valence 
quark, gluon and sea quark it is easy to see that we have n s = 2, 3 and 4 respectively. So we 
may expect (3 V ~ 3, j3 g ~ 5 and /3 sca ~ 7. 

For a rough guide to the anticipated values of the ctj parameters, we might appeal to Regge 
behaviour, since the limit x = Q 2 /2p ■ q — > corresponds to s 7P ~ 2p ■ q — > 00. In this limit the 
7P cross section is approximately proportional to 

£ elxhix) ~ (rpa^W- 1 + r^ )" 1 ) ~ (rpx 1 -^^ + rflX 1 —®). (82) 

The naive expectations are that the Pomeron and the leading secondary Reggeons have trajec- 
tories with intercepts ap(0) ~ 1.08 and cxr(0) — 0.5. The Pomeron corresponds to flavourless 
exchange so we expect the parameters a se a, 9 ~ —0.08, whereas the valence density corresponds 
to flavour exchange with a v ~ 0.5. So, in summary, we might naively expect 

xf v ~ x°-\l - x) 3 , xf 9 ~ x- a08 (l - x)\ xf sea ~ x- a08 (l - x) 7 (83) 

types of behaviour. 

In practice, the heavy quark densities, c, b, require special treatment. These are particularly 
important at small x, especially as Q 2 increases. We can see the problem by noting that for 
Q 2 ~ ml the charm quark does not act like a parton, but instead is created in the final state 
by photon-gluon fusion, 75 — > cc. On the other hand for Q 2 ^> m^, clearly c behaves like 
a massless parton. It is therefore necessary to use a variable flavour number scheme [12] in 
which we match a 3- to a 4-flavour parton description as we evolve up through the charm quark 
threshold, Q 2 ~ m 2 c . 

Table 1 highlights some processes used in the global fits, and their primary sensitivity to 
the parton densities. The kinematic ranges of the fixed-target and collider experiments are 
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Figure 14: Parton densities, xfi(x,fi 2 ), at /i 2 = 20 and 10 4 GeV 2 , obtained in a recent NNLO 
global analysis [11]. The dominance of the gluon at small x and of the valence quarks at large 
x is clearly evident. The uncertainties shown only reflect the errors of the experimental data. 
A discussion of the theoretical errors can be found in [13]. 




M (GeV) 

Figure 15: The uncertainty in the qq and gg parton luminosities for producing a state of mass M 
at the LHC, arising from the experimental errors of the data fitted in a global parton analysis. 
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Table 1: Lepton- nucleoli and related hard-scattering processes (whose data are used in the 
global parton analyses) and their primary sensitivity to the parton distributions that are probed. 





Main 


PDFs 


Process 


Subprocess 


probed 



£ ± N^£ ± X 7 *g^g g(x < 0.01), q, q 

£ + {£~)N -> v{y)X W*q -> q' 

u{u)N -> £~{£ + )X W*q->q' 

v N -> n + [i~X W*s -> c -> /i + s 

£N -> £QX 7*g -> Q Q = c,b 

i*g - < o.oi) 

-> 7X qg^iq 9 

pN — > fi + fi~X qq — > 7* g 

pp, pn — > n + fi~X uu, dd — > 7* m — d 

■urf, rfw — * 7* 

ep, en — > e7rX 7*g — ► g 

/^-^W-^^X u,d,u/d 

pp^jet+X gg,qg,qq^2j q, g(0M < x < 0.5) 



complementary (as is shown in Fig. 4) which enables the parton densities to be determined 
over a wide range in x and Q 2 . The analyses can now be done to NNLO. An example 11 of the 
resulting parton distributions is shown in Fig. 14. 

The gluon density is the most poorly known parton distribution. At small x ( < 0.01) it is 
constrained by the HERA DIS scaling violations, and for values of x up to about 0.5 by the 
Tevatron jet data. The momentum sum rule also gives an important constraint. 

Thanks to the HERA experiments, the parton densities are well-known 12 down to about x ~ 
10~ 3 . Also they are well-known up to x ~ 0.5. What are the implications of the uncertainties 13 
in the parton densities for the LHC experiments? Some idea can be obtained from Fig. 15, 
which shows the uncertainties in the C q q and C gg parton luminosities relevant to the production 
of a state of mass M at the LHC. The parton luminosities are defined as 

Cab = C ab / -±f a (x a )f b (T/x a ) (84) 

Jt x a 

11 Comprehensive sets of parton densities available as programme-callable functions can be found in 
http://durpdg.dur.ac.uk/HEPDATA/PDF. 

12 We discuss possible corrections arising from the resummation of logl/x terms and from absorptive effects, 
both of which lie outside pure DGLAP, in Sections 10 and 11 respectively. We shall see that, at low scales, the 
parton densities have large uncertainties for x < 1CP 3 . 

13 Detailed discussions of the uncertainties arising in the global analyses can be found in [14, 15, 16, 13]. 
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Figure 16: Partonic x, Q 2 domains sampled by the LHC and HERA, as well as fixed-target DIS 
experiments. The rapidity interval for the production of a Higgs boson of mass 120 GeV at the 
LHC is indicated by an open arrow; the relevant parton distributions should be reliable from 
DGLAP evolution of global analyses of HERA, fixed-target DIS, and Tevatron jet data. The 
possibility of the LHC experiments probing the region x < 1(T 4 , Q 2 > 10 GeV 2 is mentioned 
at the end of Section 12. 



where C a b is a colour factor. Since x a XbS ~ M 2 we see Xf, = rjx a where r = M 2 /s. Due to the 
factorization theorem, the cross section for the production of the state of mass M is 

a = ^ Cab &(ab -> M;s = rs). (85) 

The widening of the gg — > M error band in Fig. 15 for M > 1 TeV is due to the lack of 
knowledge of the gluon at high x. This plot does not include the theoretical errors in a pure 
DGLAP parton analysis. Nevertheless, for the predictions of the cross sections of the central 
production of high mass systems at the LHC, the uncertainty coming from parton densities is 
less than ±10%. This is also clear from an inspection of Fig. 16. 

As an example, we show in Fig. 17 the predicted cross sections[17] for production at 
the LHC. At zeroth order we only have the gg-driven subprocesses ud — > W + and du -—>■ W~; 
so we expect the parton luminosity errors to be relatively small. The cross section inequality 
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Figure 17: LO, NLO and NNLO predictions for the rapidity distribution of W production 
at the LHC The width of the bands reflects the uncertainty coming from the variation of the 
scale in the interval Miy/2 < /i < 2Myy- The NNLO prediction is the very narrow band lying 
within the NLO error band. 



a{W + ) > cr(W~) reflects u v > d v . Also note the rapid decrease in the uncertainty due to scale 
changes as we proceed from LO^NLO— >NNLO. Allowing for uncertainties from all sources, 
the W ± production cross section is predicted to an accuracy of ±5%, which enables it to be 
considered as a luminosity monitor for the LHC. 

In Fig. 18 we show the cross sections in nb for various processes at the Tevatron and at the 
LHC. If the collider luminosities were 10 33 cm _2 s _1 , then the scale on the plot also gives the 
number of events which would occur each second. Note that eventually the LHC is planned to 
achieve a luminosity some 10 times greater than this. 



10 Beyond DGLAP: low x partons and BFKL 

Fig. 19 shows the physical phenomena we expect to be appropriate in various regions of the 
log(l/x) - \og(Q 2 ) plane. We shall discuss them here and in the next Section. Overlaid is a line 
indicating the reach achieved by the HERA experiments. Of course the position of this line is 
well known, see Fig. 4. However the positions of the various domains relative to this line are 
not well established. Certainly HERA has opened up the small x domain, with DIS structure 
function measurements reaching down to x ~ 10 -4 while Q 2 is still in the perturbative domain. 

So far our approach has been to work with DGLAP evolution truncated at a fixed per- 
turbative order. This pure DGLAP approach has been phenomenologically successful, even, 
surprisingly, down to x ~ 1CT 4 with Q 2 ~ 2 GeV 2 . Nevertheless, although the global parton 
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Figure 18: The cross sections (in nb) for various processes at the Tevatron and the LHC. For 
the LHC luminosity quoted, the scale also corresponds to the number of events/second. We 
also give an indication of the physics which may be probed by the processes at the LHC. Note 
that the rates of Higgs and SUSY particle production do not include the dilution of a possible 
signal due to the branching fraction of the particular channel investigated. Moreover note 
how important it is to reduce the huge background and to overcome "pile-up" from multiple 
events per bunch crossing at the higher luminosity. Of course it would be even more exciting 
to discover something totally unexpected. 
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Figure 19: Schematic sketch of the physical phenomena in various regions of the log(l/x) - 
log(Q 2 ) plane, compared to the kinematic reach of HERA. The gluonic content of the proton, 
as resolved by a Q 2 probe, is also indicated. DGLAP evolution takes us up in Q 2 and so the 
partonic constituents are resolved more finely. The BFKL equation takes us to small x, with 
the gluon density xg growing as x~ x , but the resolution in the transverse plane remaining at 
approximately l/Q. As x decreases, the partonic content increases, and at some stage the 
partons recombine (absorptive effects), and eventually saturate. 

analyses describe the data satisfactorily in this regime 14 , it does not mean that the parton 
distributions are reliable here. We know pure DGLAP is incomplete at small enough x. 

To explore the small x regime, we first note that DGLAP is equivalent to assuming that 
the dominant dynamical mechanism leading to DIS scaling violations is the evolution of parton 
emissions strongly-ordered in transverse momenta. However, at small x the evolution occurs 
over large rapidity intervals (~ lnl/x). The higher-order corrections to the splitting (and 
coefficient) functions contain one additional power of lnl/x for each additional power of a s . If 
we keep just the leading lnl/x terms then the small x behaviour of the P gg splitting function, 
for example, has the form 

xP gg {x) -> Alo a.s + A 2 i afilal/x + A 32 a 3 s ln 2 l/x + A A3 a 4 ln 3 l/x + ... , (86) 

14 The gluon has a valence-like behaviour, although the uncertainties are large in this domain. Nevertheless 
its behaviour is quite different to the growth of the sea-quark distributions as x — > 0. Such a result looks strange 
from the Regge viewpoint where the same vacuum singularity (Pomeron) should drive both the sea quarks and 
the gluons; i.e. the same small x behaviour is expected for sea quark and gluon distributions. 
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whereas in NNLO DGLAP, for example, P gg contains only the terms up to a 3 s . Clearly, at 
small x, when a s Inl/x ~ 1, a resummation of all of the terms in the series is necessary. 
The resummation of the leading log (LL X ) terms, a™log n_1 l/:r, is accomplished by the BFKL 
equation 15 . The BFKL equation [18] will be discussed in detail in the lectures of Victor Fadin 
[19], Lev Lipatov [20] and Al Mueller [21]. Here we will just include some introductory remarks. 

At low x we have diffusion or "random walk" in the logarithm of transverse momenta as we 
proceed along the emission chain. We no longer have the strong ordering in kt which is true 
in DGLAP evolution. For this reason the BFKL equation is for the gluon density, f(x,k^), 
unintegrated over k t . Recall that the gluon dominates at low x. The BFKL equation has the 
structure 

df 

gtffcj = **/. = V. (87) 
which, at small x, has the the solution 

/, ~ e^V*) „ x -x „ (88) 

where A = 12o; s ln2/7r is the leading eigenvalue of the BFKL kernel K. This has an analogous 
form to the Regge-pole exchange behaviour of the amplitude, 

A(s,t) ~ £/?*(*)(-) > (89) 

R W 

which is the cornerstone of the description of high-energy "soft" hadron-hadron interactions; 
dji(t) is the trajectory of Reggeon R in the complex angular momentum plane. For colour- 
octet exchange the BFKL equation describes a Reggeized gluon with trajectory oi g (t), while for 
colour-singlet exchange, which is relevant to this discussion, it leads to a cut in the complex 
angular momentum, j, plane corresponding to two Reggeized gluons exchanged - often called 
the perturbative Pomeron. Note that the generalized gluon distribution f g corresponds to the 
two-gluon exchange amplitude. Its behaviour at low x is driven by the rightmost singularity 
(branch point), j = 1 + A, produced by the two-gluon cut, where the value of A obtained from the 
BFKL equation is given above. Since the behaviour of f g is driven by a cut (and not an isolated 
pole) in the j-plane, a prefactor l/\/\ns will appear in (88). The possible connection between 
(89) and (88) is indicated by a horizontal 'block' arrow in Fig. 19. In the "soft" regime the 
hadrons are Reggeized, while in the perturbative QCD BFKL regime the constituent partons 
are Reggeized. How to go from one regime to the other has not been solved. For example, 



15 DGLAP and BFKL arc different limits of a more general evolution of parton densities, which is an ordered 
evolution in the angles of the emitted partons. At LO we have strong ordering of the emission angles, <C 
on the other hand if, at one step of the evolution 0i <~ 0i+i, then this contribution is included inside 
the NLO splitting function. In the collinear approximation of DGLAP the angle increases due to the growth 
of the transverse momentum kt, while in BFKL the angle (9 ~ kt/k\\) grows due to the decreasing longitudinal 
momentum fraction as we proceed along the emission chain from the proton. Introductory discussions of the 
BFKL equation can be found, for example, in Refs.[25, 26]. 
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what is the relation 16 of the BFKL or perturbative QCD 'Pomeron' (given by a ladder diagram 
formed from the exchange of two t-channel Reggeized gluons) to the 'Pomeron' describing soft 
high-energy proton-proton interactions? 

Coming back to the discussion centred on the perturbative expansion of equation (86), 
we note that in the small x region the gluon dominates, and only P gg and P qg contain LL X 
contributions. These are positive but smaller than naively expected; it turns out that A 2 \ = 
A 32 = 0, and even A 54 = 0, in (86). Now the next-to-leading (NLL X ) terms, a™\og n ~ 2 l/x, have 
also been calculated [23]. These give a large, negative, contribution to the gluon, leading to 
instability at small x. In fact A of (88) is now given by 

A = 12a s ln2/7r (1 - 6.5a s ). (90) 

This problem has been the subject of considerable investigation. Clearly, the higher-order 
contributions, NNLL.J,, NNNLL X ,... are important. However it took about 10 years to calculate 
the NLLa, contributions, so to compute the next order or two appears unrealistic, and even 
then may not converge to a stable result. Instead, the procedure that has been followed is to 
identify a few physical QCD effects that lead to large higher-order corrections and then to resum 
them. Indeed, this all-order resummation of the main effects is found to tame the wild (LL^ 
— > NLLa;) behaviour; a readable review is given in [24]. The approaches of the various groups 
have reached similar conclusions: the approximate all-order resummed BFKL framework leads 
to the behaviour that 

xg ~ a; -0 - 3 as x -> (91) 

at low scales 17 . 

In practice, it is found that this power-like growth only sets in at very small x. In terms of 
DGLAP evolution all the BFKL effects should be included in the resummed splitting functions 
used to describe the transition between two quite different scales, that is between partons whose 
transverse momentum are very different. In such a case the power growth (91) will be included 
in the resummed P gg . However the resummed xP gg has a dip centred at x ~ 10~ 3 , and the 
power growth is only evident below x ~ 10~ 5 . Indeed the resummed xP gg and the NNLO 
DGLAP xP gg are in good agreement down to x ~ 10~ 3 . 

To make quantitative predictions in the small x domain, x < 10~ 4 with Q 2 ~ 2 GeV 2 , where 
no data exist, is extremely difficult. We need to obtain the resummed ln(l/x) solution starting 
from some non-perturbative amplitude at Q = Qq. This non-perturbative distribution (which 
is analogous to the 'input' in the DGLAP approach) is not known theoretically. Either one has 
to fit it to data (but again low x data are needed) or to use some phenomenological model (for 
example, based on a Regge parametrization). 

16 See [22] for a phenomenological study. 

17 What do the data say? If the F 2 data are fitted to the form x~ x for x < 0.01, then it is found that A 
grows approximately linearly with logQ 2 from A ~ 0.1 passing through A = 0.3 at Q 2 ~ 40 GeV 2 . The simple 
assumption that this reflects the behaviour of the gluon, with F2 driven by the g — > qq transition is much too 
naive. Indeed the global analyses give a gluon which is valence-like at small x at the input scale. 
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We conclude that the parton densities are unknown in the region x < 1CT 4 . At very small 
x we have the estimate that gluon density might behave as xg ~ x~ x with A ~ 0.3. However, 
as x decreases, at some stage this behaviour will violate unitarity. Here the recombination of 
gluons (absorptive effects) come to the rescue, and tame the violations of unitarity. To this we 
now turn. 



11 Absorptive effects 

The saturation of parton densities (A = 0) may be obtained using the Gribov-Levin-Ryskin 
(GLR) equation [27] or the more precise Balitski-Kovchegov (BK) equation [28]. These equa- 
tions sum up the set of so-called fan diagrams which describe the rescattering of intermediate 
partons on the target nucleon. The screening caused by these rescatterings prohibits the power 
growth of the parton densities. 

The GLR equation for the gluon may be written in the symbolic form 

The non-linear shadowing term, — [g] 2 , describes the recombination of gluons. It arises from 
perturbative QCD diagrams which couple Ag to 2g — that is two gluon ladders recombining 
into a single gluon ladder, which is called a fan diagram. The minus sign occurs because 
the scattering amplitude corresponding to a gluon ladder is predominantly imaginary. The 
parameter R is a measure of the transverse area nR 2 where the gluons are concentrated. 

The BK equation is an improved version of the GLR equation. It accounts for a more precise 
form of the triple-pomeron vertex and can be used for the non-forward amplitude. The GLR 
equation, based on DGLAP evolution, was in momentum space; whereas the BK equation, 
based on the BFKL equation, is written in coordinate space in terms of the dipole scattering 
amplitude iV(x, y,Y) = N xy (Y). Here x and y are the transverse coordinates of the two 
t-channel gluons which form the colour-singlet dipole, and Y = ln(l/x) is the rapidity. The 
BK equation reads 

dN*y = 3a 1 rd 2 z (x-y) 2 

dY n J 2tt (x - z)2(y - z)2 iiVxz + iVyz iVxy iVxz iVyz| ' ^ 

where, interestingly, the non-linear and linear terms have the same BFKL kernel K, which is 
shown explicitly in (93). For small dipole densities, N, the quadratic term in the brackets may 
be neglected, and, indeed, (93) reproduces the conventional BFKL equation. However for large 
N, that is iV — > 1, the r.h.s. of (93) vanishes, and we reach saturation when N — 1. The 
equation sums up the set of fan diagrams where at the lower (small Y) end the target emits 
any number of pomerons (i.e. linear BFKL amplitudes), while at the upper (large Y) end we 
have only one BFKL dipole. 
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In principle, it would appear more appropriate to use the BFKL-based BK equation to 
describe the parton densities at low x. It is an attempt to describe saturation phenomena. 
However it is just a model and cannot, at present, be used to reliably estimate absorptive 
effects at small x. 

Is there any evidence of the onset of absorptive effects in the experimental data? These 
should occur first at low x and low Q 2 , see Fig. 19. However, there is no conclusive evidence that 
absorptive effects are important in the HERA data in the perturbative regime, Q 2 > 1 GeV 2 . 
The various claims that are frequently made have been recently comprehensively discussed in 
[29]. It is seen that none of them, including the observed 'flat' ratio of (diffractive DIS/inclusive 
DIS) or the observation of geometric scaling, provide any compelling evidence of saturation 
effects. 

Of course, as x decreases we know that ultimately absorptive effects must be present. In 
principle, we should be able to estimate their contribution from knowledge of the structure 
functions for diffractive DIS, via 

AF 2 abs ~ -F 2 D , (94) 

where F 2 D is the structure function for the process 7*p — > X +p in which the slightly deflected 
proton and the cluster X of outgoing hadrons are well separated in rapidity [30]. 

12 Conclusions 

The great improvement in the precision and range of deep inelastic and related hard scat- 
tering data over the last few years has enabled the partonic structure of the proton to be 
well determined in the 1CT 3 < x < 0.5 interval, so we are able to make reliable predictions for 
the production of new massive states at the LHC. Global analyses are now available at NNLO. 
These analyses require particularly careful treatment at the heavy quark thresholds, see Thorne 
[12] and references therein. A surprise is that a pure DGLAP description is able to describe 
all features of the data down to Q 2 = 2 GeV 2 , in spite of the fact that the global fits are now 
quite tightly constrained. The allowance of beyond-DGLAP effects is not found to improve the 
description, see, for example [13]. 

Another surprise is that the global analyses reveal that the gluon has a valence-like small 
x behaviour at the low input scale, unlike the sea quark distribution which behaves as ex- 
pected, see (83). Such a result looks strange from the Regge viewpoint where the same vacuum 
singularity (Pomeron) should drive both the sea quarks and the gluons; i.e. the same power 
A 9 = A sea is expected for sea quarks and gluons. Note that global analyses are only reliable 
down to x ~ 10~ 3 . 

The low x domain, x < 10~ 4 , is unchartered territory. Is it possible for the LHC experiments 
to determine the behaviour of partons in the important low x region below 10~ 4 at low scales? 
One possibility is Drell-Yan production in which events are observed with the n + ji~ 
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invariant mass as low as possible and the rapidity as large as possible. For example, for 
Mjjj, = 4 GeV and y m = 3, we sample quarks at x — 1.4 x 10~ 5 , see Fig. 16. This process 
samples predominantly the sea quark distributions. To study the small x behaviour of the 
gluon at low scales we may consider Xc production 18 , or prompt photon production driven by 
the subprocess gq — > 79, or perhaps gg — > bb. These studies may also require an improvement 
in theoretical formalism. 
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